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FLUCTUATIONS  OF  THERMODYNAMIC  PROPERTIES  OF  SUPERCOOLED  LIQUID  WATER 

G.  E.  Walrafen  and  W.-H.  Yang* 

Chemistry  Department,  Howard  University,  Washington,  DC  20059 


It  has  been  known  for  many  years  that  fluctuations  of  several  thermo¬ 
dynamic  properties  of  supercooled  liquid  water  are  highly  unusual. 

An  example  is  the  heat  capacity  at  constant  pressure,  Cp.  As  water 
is  supercooled,  Cp  rises  very  rapidly,  particularly  below  about 
-20°C  ( Angel 1  1982).  Because  S2-(S)2  =  kBCp,  where  kB  =  1.3806  x 10"23 
J°K"1,  it  follows  that  the  mean  square  fluctuations  of  the  entropy, 

S,  must  .increase  as  the  temperature  decreases .  This  situation  is 
contrary  to  intuition,  and  immediately  signals  the  presence  of 
unusual  effects,  e.g.,  it  suggests  a  continuous  transformation 
between  structures  of  markedly  different  entropy  (Speedy  1984). 
Fluctuations  of  the  enthalpy,  energy,  and  molal  volume,  which  increase 
with  decreasing  temperature__  at  constant  pressure  (T,p,N  ensemble) 
may  be  obtained  from:  H2-(H)2  =  kBT2Cp  =  kBT2 [Cv+( b^TV/ioj) ] , 

E2_(E)2  =  kBT{TCv  +  [OE/3V)t12<tv>»  and  V2-(V)2  =  kBTiexV  (Kell  1982; 
Kubo  1981;  McQuarrie  1976).  (Temperature  fluctuations,  T2-(T)2  = 
kBT^fcv,  increase  at  an  increasing  rate  between  -30°C  and  150°C, 
and  pressure  fluctuations,  p2-(p)2  =  kBT/Vic5,  increase  to  a  maximum 
near  100°C,  and  then  decline  with  increasing  temperature  to  150°C.) 

B  =  V-M3V/3T)  ,  kt  -  -V-l(3V/3p)T,  <s  =  -V”1  ( 3V/3p )  -  =  icxCv/CD,  and 
(3E/3V)t  =  Pi  =  (BT/kt)-P.  V 

Fluctuations  corresponding  to  the  canonical  or  T,V,N  ensemble,  for 
which  V  is  fixed,  are  also  of  interest  for  liquid  water.  For  this 
ensemble,  E2-(E) 2  =  kBT2cv  and  (E-Ep  =  k2{T4  (  3CV/3T) v  +  lT3Cy) 

(Kubo  1981).  E2-(£j2  (T,V,N)_is  important  because  it  constitutes 

the  leading  term  in  both  H2-(H)2  and  E2-(e)2  (T,p,N).  Hence,  _ 

E2- ( E ) 2  ^T,V,N)  may  be  regarded  as  a  "baseline"  above  which  H2-(H)2 
and  E2- ( E ) 2  (T,p,N)  (which  are  equal  for  practical  purposes)  both 
rise,  below  and  above  the  temperature  of  minimum  molal  volume,  4°C. 

A  plot  of  E2-(E)2  (t,V,N)  versus  temperature  from  -30°C  to  150°C  was 
found  to  be  so  nearly  perfectly  linear  that  small  negative  deviations 
near  -30°C  suggested  error  in  Cy  (determined  from  Cp  -  B2TV/VX ) , 
rather  than  a  failure  of  the  linearity.  Therefore,  if  E2-(E)2  = 
a'  +  b'T,  it  follows  that  CVT2  =  a  +  bT,  where  b>0  and  a<0.  Cy 
values  calculated  between  0°  and  100°C  were  considered  to  be  more 
accurate  than  those  over  the  remaining  parts  of  the  -30  to  150°C 
region.  Least  squares  treatment  (0°  to  100°C)  yielded  the  empirical 
relation,  CV  =  -4 , 691 , 495T“2  ♦  3 7 , 930 . 92T"1 ,  in  J°k-lmol-1.  (Cp 
and  Cy  fcannot  be  negative,  but  this  is  irrelevant  for  terms  in  an 
empirical  equation.)  To  test  the  empirical  equation  values  of 
Ex-E277.16  =  EfE4  were  obtained  from  Cy  =  Cp  -  (B2TV/kx)  by  graphi-  ,, 
cal  integration  versus  the  Cy  values  from  the  equation,  obtained  by 
analytical  integration.  The  comparison  is  shown  below.  Seme  large 
deviations  are  seen  to  be  0.25%  at  -20°C,  where  Cp  is  given  only  to 
3  significant  figures  (Angell  1982);  and  0.1%  at  -30°C,  where  B  is 
given  only  to  2  significant  figures  (Kell  1982). 
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The  third  moment,  (E-E)3,  is  readily  calculated  from  Cy  =  aT-2  + 
bT”^»  namely  ( E-E ) J  =  bksT.  This  moment  is  a  measure  of  the 
skewness,  which  is  relatively  very  small,  as  seen  from 
(E-eH/E2-(E)2  =  bkB/Cv  =  6.92  x  10-21,  at  4°C. 

Speedy  (1984  )  has  modelled  the  fluctuation  aH  aV  in  terms  of 
pentagonal  rings  of  hydrogen-bonded  water  molecules  thought  to 
exist  in  supercooled  water.  These  rings  are  consistent  with  the 
presence  of  bulky  polyhedral  species.  Speedy  (1984)  formulated 
aHaV  as  keST^V,  which  is  the  leading  term  in  the  formulation  of 
Kubo  (1981)  A EAV  =  kBT{T3V-pV<T } .  (H  is  nearly  the  same  as  E  for  __ 
liquid  water.)  In  view  of  the  present  evidence  that  S 2-(S)2,  H2-(H)2 
E2-(E)2  and  V2-(V)2  all  increase  as  the  temperature  of  supercooled 
water  decreases ,  it  seems  nearly  certain  that  special,  bulky,  fully- 
hydrogen-bonded  structures  exist  in  highly  supercooled  water,  where 
the  molal  volume  rises  sharply.  Thus,  these  structures  would  be 
volume  modulated ,  whereas  the  breakage  of  hydrogen  bonds  is 
temperature  modulated  (Walrafen  1986).  Hence,  these  two  modulation 
mechanisms  must  be  considered  for  water. 

Finally,  one  simple  calculation  highlights  the  great  difference 
between  supercooled  water  and  ice  at  the  same  temperature.  For 
supercooled  water  at  -13°C,  (3E/3V)T  =  py  =  -1900  bar,  whereas  the 

corresponding  py  for  ice  ( 5>0  for  ice,  8<0  for  water,  -13°C)  is 
about  +3000  bar  (Leadbetter  1965).  This  qualitative,  large  disparity 
of  nearly  5000  bar  emphasizes  the  differences  between  ice  and  super¬ 
cooled  water.  Large  bulky  polyhedral  structures,  e.gr.,  clathrate- 
type  structures,  do  not  exist  in  ice  Ih,  but  such  structures  may 
exist  in  supercooled  liquid  water. 
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